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Abstract 

■ Based on the reduced density matrix method, we compare two dif- 

Q\ ' ferent approaches to calculate the dynamics of the electron transfer in 

ON ■ systems with donor, bridge, and acceptor. In the first approach a vibra- 

tional substructure is taken into account for each electronic state and the 
O ■ corresponding states are displaced along a common reaction coordinate. 

In the second approach it is assumed that vibrational relaxation is much 
faster than the electron transfer and therefore the states are modeled by 
' electronic levels only. In both approaches the system is coupled to a bath 

of harmonic oscillators but the way of relaxation is quite different. The 
theory is applied to the electron transfer in H 2 P — ZnP — Q with free-base 
porphyrin (EbP) being the donor, zinc porphyrin (ZnP) being the bridge 
qq | and quinone (Q) the acceptor. The parameters are chosen as similar as 

\Q . possible for both approaches and the quality of the agreement is discussed. 
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1 Introduction 



Long-range electron transfer (ET) is a very actively studied area in chemistry, 
biology, and physics; both in biological and synthetic systems. Of special in- 
terest are systems with a bridging molecule between donor and acceptor. For 
example the primary step of charge separation in the bacterial photosynthe- 
sis takes place in such a system JlJ. But such systems are also interesting for 
q-i synthesizing molecular wires It is known that the electronic structure of 

the bridge component in donor-bridge-acceptor systems plays a critical role 
g ||. When the bridge energy is much higher than the donor and acceptor 
energies, the bridge population is close to zero for all times and the bridge site 
just mediates the coupling between donor and acceptor. This mechanism is 
called superexchange and was originally proposed by Kramers pi to describe 
the exchange interaction between two paramagnetic atoms spatially separated 
by a nonmagnetic atom. In the opposite limit when donor and acceptor as well 
as bridge energies are closer than ~ k&T, the bridge site is actually populated 
and the transfer is called sequential. The interplay between these two types of 
transfer has been investigated theoretically in various publications g, [7|, §, § . 

In the present work we compare two different approaches based on the re- 
duced density matrix formalism. In the first model one pays attention to the 
fact that experiments in systems similar to the one discussed here show vibra- 
tional coherence [10, 11]. Therefore a vibrational substructure is introduced for 



each electronic level within a multi- level Redfield theory [12, 13]. Below we call 
this the vibronic model. In the second approach only electronic states are taken 
into account because it is assumed that the vibrational relaxation is much faster 
than the ET. This model is referred to as tight-binding model below. In this 
case only the relaxation between the electronic states remains. Such a kind of 



relaxation has been phenomenologically introduced for ET by Davis et al. [14] 
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and very recently derived in our group |l5|, [l6| as a second order perturbation 
theory in the system-bath interaction similar to Redfield theory. The vibronic 
and the tight-binding model are described in the next section and compared in 
Section 3. 



2 Theory 

For the description of charge transfer and other dynamical processes in the 
system we introduce the Hamiltonian 

H = H S + H B + H SB , (1) 

where Hs denotes the relevant system, Hb the dissipative bath, and Hsb the 
interaction between the two. Before discussing the system part of the Hamil- 
tonian in Sections 2.1 and 2.2, we describe the bath and the procedure how to 
obtain the equations of motion for the reduced density matrix, because this is 
the same for both models studied below. The bath is modeled by a distribution 
of harmonic oscillators and characterized by its spectral density J{u). Starting 
with a density matrix of the full system, the reduced density matrix of the 
relevant (sub)system is obtained by tracing out the bath degrees of freedom 
While doing so a second-order perturbation expansion in the system-bath 
coupling and the Markov approximation have been applied JT^]. 



2.1 Vibronic model 

The bridge ET system H2P — ZnP — Q with free-base porphyrin (H2P) being 
the donor, zinc porphyrin (ZnP) the bridge, and quinone (Q) the acceptor is 
modeled by three diabatic electronic potentials, corresponding to the neutral 
excited electronic state |1) = IH2P* — ZnP — Q), and states with charge sepa- 
ration |2) = |H 2 P + - ZnP" - Q), |3) = |H 2 P + - ZnP - Q~) (see Fig. 1). Each 
of these electronic potentials has a vibrational substructure. The vibrational 
frequency is assumed to be 1500 cm -1 as a typical frequency within carbon 
structures. The potentials are displaced along a common reaction coordinate 
which represents the solvent polarization |l8|]. Following the reasoning of Mar- 
cus jyj the free energy differences AG mn corresponding to the electron transfer 
from molecular block n to m (n = 1, m = 2, 3) are estimated to be [19, 2C] 

red e e 2 1 

AG mn = — — E ex — h AG mn (e s ) (2) 

47re e s r mn 

with the term AG mn (e s ) correcting for the fact that the redox energies E™ and 
-E^ ed are measured in a reference solvent with dielectric constant e r s ei : 

e 2 ( 1 1 \ / 1 1 



The excitation energy of the donor H2P — ► H2P* is denoted by E ex . r n denotes 
the radius of either donor (1), bridge (2), or acceptor (3) and r mn the distance 
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between two of them. They are estimated to be t\ = T2 = 5.5 A, r% = 3.2 A, 
ris = 12.5 A, and r 13 = 14.4 A ||, |0|. 

Also sketched in Fig. 1 are the reorganization energies X mn = X\ nn + A^ n . 
These consist of an internal reorganization energy A^ n , which is estimated to 
be 0.3 eV |20[| , and a solvent reorganization energy |]l^] 

^„ = ^( J_ + _L __L)(±_I). (4) 

47re \2r m 2r n r mn / Ve^ e s / 

Further parameters are the electronic couplings between the potentials. First 
it should be underlined that V\% = because of the spatial separation of H2P 
and Q. So there is no direct transfer between donor and acceptor. The other 
couplings are V12 = 65 meV and V23 = 2.2 meV pQ[| . The damping is described 
by the spectral density J(u) of the bath. This is only needed at the frequency 
of the vibrational transition and is determined J(w v ib)/Wib = 0.372 by fitting 
the ET rate for the solvent methyltetrahydrofuran (MTHF). In the vibronic 
model the spectral density is taken as a constant with respect to e s . 

Next the calculation of the dynamics is sketched. Starting from the Liou- 
ville equation, performing the abovementioned approximations the equation of 



motion for the reduced density matrix p^ u can be obtained [12 



The index /j, combines the electronic quantum number m and the vibrational 
quantum number M of the diabatic levels E^. v^ v = V mn Fpc('m', M, n, N) com- 
prises Franck-Condon factors Fpc and the electronic matrix elements V mn . The 
third term describes the interaction between the relevant system and the heat 
bath. Equation (||) is solved numerically with the initial condition that only 
the donor state is occupied in the beginning. The population of the acceptor 
state 

Ps(t) =I>3Af3Af(t) (6) 



M 



and the ET rate 



P 3 (oo) m 

fc ET = 00 ( 7 ) 

fdt(l-P 3 (t)) 




are calculated by tracing out the vibrational modes. 



2.2 Tight-binding model 

The reasoning for the following system Hamiltonian is the assumption that the 
vibrational excitations are relaxed on a much shorter time scale than the ET 
time scale. Therefore only electronic states without any vibrational substruc- 
ture are taken into account (see Fig. 2). As a consequence the relaxation during 
the ET process has to be described in a different manner than in the previous 
subsection. If now relaxation takes place, it takes place between the electronic 
states and not between vibrational states within one electronic state potential 
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surface. A similar model has been introduced phenomenologically by Davis et 
al. jL4|] who solved it in the steady state limit. 

The energies of the electronic states E m are chosen to be the ground states 
of the harmonic potentials given in the previous section. So they vary with the 
dielectric constant. The electronic coupling is fixed in two different ways. In the 
naive way they are chosen to be the same as in the vibronic model. But because 
in the tight-binding model there is no reaction coordinate, in a second version we 
scale the electronic couplings with the Franck-Condon overlap elements between 
the vibrational ground states of each pair of electronic surfaces 

Vmn = V m nFFc(m, 0, n, 0) = V mn exp j^ mn l . (8) 

In the vibronic model not only the free energy differences AG but also the 
reorganization energies A scale with the dielectric constant e s . Due to this 
scaling of A the system-bath interaction is scaled with the dielectric constant 



e s . In the high temperature limit the reorganization energy is given by [21] 



A^r^. ( 9) 

Jo oj 

This relation is taken as motivation to scale the tight-binding spectral density 
with e s like the reorganization energies A in the vibronic model. In the present 
calculations I^i = 1^3 = T is assumed. The absolute value of the damping rate 
r between the electronic states (see Fig. 2) is then determined by fitting the 
ET rate for the solvent MTHF to be T = 2.8 x 10 11 s -1 . 

The advantage of the tight-binding model is the possibility to determine the 
transfer rate &et an d the final population of the acceptor state either numeri- 
cally or analytically. We employ the rotating wave approximation because we 
are only interested in the reaction rates here. For the analytic calculation three 
extra assumptions have to be made: small bridge population, the kinetic limit 
t 3> r _1 , and the absence of initial coherence in the system. But for all situ- 
ations described in this paper the differences between analytic and numerical 
results without the extra assumptions are negligible. The analytic expressions 
are 

, . g23(g!2 ~ 932) , 1fV . 
ket = 923 H ; (10) 

521 + 923 

and 

D / \ 912923 /, x-1 ni N 

P 3 (oo) = ■ (Atet) , (11) 

921 + 923 

which contain both, dissipative and coherent contributions 

J2(dmk + dkn) 

9mn — d mn H p - " (1^) 



^ 2 \ ^mn + \ 



J2(dmk + dkn) 



Herein the d m n are just abbreviations for r mn |n(o;r nn )| and n(uj m n) denotes the 
Bose distribution at frequency u> mn = (E m — E n )/h. For details and comparison 



with the Grover-Silbey theory [22] as well as the Haken-Strobl-Reineker theory 



1 23] we refer the reader to Ref. |ll 
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3 Comparison 



In Fig. 3 it is shown how the minima of the potential curves change with varying 
the solvent due to the changes in Eqs. (Q) to (|J). The solvents are listed in 
Table 1 together with their parameters and the results for the ET rates in 
both models. For larger e s the coordinates of the potential minima of bridge 
and acceptor increase while their energies decrease with respect to the energy 
of the donor. The energy difference between donor and bridge decreases with 
increasing e s . This makes a charge transfer more probable. For small e s the 
acceptor state is higher in energy than the donor state; nevertheless there is a 
small ET rate due to coherent mixing. 

For fixed eoo the ET rate is plotted as a function of the dielectric constant 
e s in Fig. 4. The ET rate in the vibronic model increases strongly for small 
values of e s while the increase is very small for e s in the range between 5 and 
8. The increase for small values of e s is due to the fact that with increasing e s 
the minimum of the acceptor potential moves from a position higher than the 
minimum of the donor level to a position lower than the donor level. So the 
transfer becomes energetically favorable. This can also be seen when looking at 
the results for the tight-binding model without scaling the electronic coupling 
with the Franck-Condon factor. In this case the ET rate increases almost 
linearly with increasing e s . The effect missing in this model is the overlap 
between the vibrational states. If one corrects the electronic coupling in the 
tight-binding model by the Franck-Condon factor of the vibrational ground 
states as described in Eq. (^), good agreement is observed between the vibronic 
and the tight-binding model. 

The ET rate for the vibronic model shows some oscillations as a function 
of e s . This is due to the small density of vibrational levels in this model with 
one reaction coordinate. All three electronic potential curves are harmonic and 
have the same frequency. So there are small maxima in the rate when two 
vibrational levels are in resonance and minima when they are far off resonance. 
Models with more reaction coordinates do not have this problem nor does the 
simple tight-binding model. If these artificial oscillations would be absent, the 
agreement between the results for the tight-binding and the vibronic model 
would be even better, because the rate for the vibronic model happens to have 
a maximum just at the reference point e s = 6.24 which we have chosen to fix 
the spectral density, i. e. for MTHF. 

The comparison of the two models has been made assuming that the scaling 
of energies as a function of the dielectric function is correct in the Marcus 
theory. There have been a lot of changes to Marcus theory proposed in the last 
years. Marcus theory assumes excess charges within cavities surrounded by a 
polarizable medium and there one only takes the leading order into account. 
Higher order terms are included in the so called reaction field theory (see for 
example p|| ) . But to compare different solvation models is out of the range of 
the present investigation. Some more details on this issue for the tight-binding 
model are given in Ref. [16]. Here we just want to note in passing that the 
effect of scaling the system-bath interaction with e s , as assumed in the present 
work for the tight-binding model, has no big effect on the ET rates. 
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As conclusion we mention that one gets good agreement for the ET rates 
of the models with and without vibrational substructure, i. e. the vibronic and 
the tight-binding model, if one scales the electronic coupling with the Franck- 
Condon overlap matrix elements between the vibrational ground states. The 
advantage of the model with electronic relaxation only is the possibility to derive 
analytic expressions for the ET rate and the final population of the acceptor 
state. But of course for a more realistic description of the ET transfer process in 
such complicated systems as discussed here, more than one reaction coordinate 
should be taken into account. Work in this direction is in progress. 
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solvent 




e s 




AG 21 


AG 3 i 


A 21 


A 31 


r 


jLel 
"'ET 


i,vib 
"'ET 










[eV] 


[eV] 


[eV] 


[eV] 


[10 11 s- 1 ] 


[10 8 s" 1 ] 


[10 8 S" 1 ] 


1. cyclohexane J20f 


2.02 


2.00 


0.976 


0.393 


0.007 


0.012 


0.042 


0.181 


0.7 


2. toluene fl25|| 




2.38 


2.24 


0.867 


0.202 


0.039 


0.069 


0.227 


1.04 


0.8 


3. anisole |26| 




4.33 


2.29 


0.590 


-0.281 


0.300 


0.524 


1.751 


4.24 


2.30 


4. dibromoethane [26 


1 


4.78 


2.37 


0.558 


-0.336 


0.312 


0.544 


1.817 


4.63 


2.45 


5. chlorobenzene 


25] 


5.29 


1.93 


0.529 


-0.388 


0.481 


0.839 


2.804 


3.21 


3.63 


6. MTHF @ 




6.24 


2.00 


0.486 


-0.462 


0.497 


0.868 


2.900 


3.59 


3.58 


7. methyl acetate 




6.68 


1.85 


0.471 


-0.489 


0.571 


0.996 


3.328 


2.96 


4.15 


8. trichloroethane 




7.25 


2.06 


0.454 


-0.512 


0.508 


0.887 


2.960 


3.98 


3.50 


9. dichloromethane |2C| 


9.08 


2.03 


0.413 


-0.590 


0.559 


0.977 


3.264 


4.00 


3.80 



Table 1: Parameters and obtained transfer rates for different solvents. The 
references behind the names of the solvents cite the sources of e s and eoo . MTHF 
stands for methyltetrahydrofuran. T denotes the damping rate in the tight- 
binding model. The ET rate for the solvent MTHF has been used to fix the 
damping parameter of the models. The reaction rates fc|J T were obtained using 
Eq. ([ji]) within the tight-binding model and the reaction rates within the 
vibronic model. 



S 




0.1 0.2 

reaction coordinate q [A] 



Figure 1: Electronic potentials and parameters of the vibronic model. 
The donor surface IH2P* — ZnP — Q) is given by the solid line, the bridge 
|H 2 P + - ZnP" - Q) by the dashed line, and the acceptor |H 2 P+ - ZnP - Q~) 
by the dotted line. 
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-0.1 0.1 0.2 

reaction coordinate q [A] 



Figure 3: Variation of the potential minima for different solvents. Squares de- 
note the bridge minima, circles the acceptor minima. The numbers correspond 
to the ordinal numbers in Table 1. The potentials are shown for solvent 6 
(MTHF). 
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23456789 
static dielectric constant e s 

Figure 4: Transfer rate as a function of the dielectric constant e s for both models 
together with experimental results pC|]. The rates for the vibronic model are 
given by the circles. The dashed line shows the rate for the tight-binding model 
with electronic couplings V mn as in the vibronic model. The solid line represents 
the rate for the tight-binding model with v mn scaled as given in Eq. (||). 
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